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Abstract 

The criticisms in manuscript arXiv:0803.4371 do not change the renormalization group equation 
for the "Viowfc" interaction nor do they affect any results in the literature. Several other potentially 
misleading statements about low-momentum interactions are also addressed. 
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A recent paper |l[ makes various statements about renormalization group (RG) equations 
for low-momentum interactions ("Vowfc") as derived by Bogner et al. in Ref. |2j. The earlier 
posted versions of the manuscript (vl and v2 of Ref. claimed that the RG equations 
were incorrect because of a missing term. After we communicated to the author an explicit 
demonstration that this term vanishes, as well as a straightforward proof that the conven- 
tional model-space methods must obey the same RG equation, the manuscript was modified 
(v3 and v4 of Ref. to criticize only the derivations in Ref. [2[ rather than the final RG 
equation or the equivalence of the model-space and RG approaches. While we agree that 
the original derivations in Ref. 0] lack formal mathematical rigor, we feel that readers of 
Ref. [l| could still be misled that there are problems with Vowfc- Therefore we have written 
this comment. 

We emphasize in particular the following point: None of the calculations in the literature 
involving Vowfc are affected by the discussion in Ref. f^j. Low- momentum nucleon-nucleon 
interactions can be derived and constructed in energy-independent form either using model- 
space methods (such as Lee-Suzuki or Okubo) or through an RG treatment. To our knowl- 
edge, there is no question about the validity of the former; the issue in Ref. jl[ concerned 
the precise equivalence of the approaches. However, in v3 and v4 the author acknowledges 
that this equivalence is indeed correct. In practice, most calculations using Viowfc have used 
the model-space methods rather than the differential equations of the RG because they are 
more robust numerically. This is true for both sharp and smooth cutoff versions of Vowfc- 
Thus doubts in Ref. [l| about the final RG equation have no impact on prior results. 

Nevertheless, the equivalence of the model-space and RG methods is important concep- 
tually. Reference proposed four derivations, none of which is rigorous in a mathematical 
sense. The most straightforward issue pointed out in Ref. [lj, which could undermine two 
of the derivations, is that the use of the completeness relation did not separately treat the 
bound states, as is usually required. We reproduce our original argument here that such a 
treatment, as carried out in Ref. does not in the end lead to changes in the RG equation 
proposed in Ref. |2j. We would welcome a more rigorous mathematical treatment of this 
and other aspects of the Viowfc RG equations, but leave that to others. 

The low-momentum partial- wave interaction Viowfc [3] is defined by the following equation 
for the half-on-shell (HOS) T-matrix (in units where h 2 /m = 1), 

Ti owk (k,k;k ) = Vi owk (k,k) + - / p dp — , (1) 

vr Jo k p 2 

where Cauchy principal value integrals are implied throughout and the quantum numbers 
labeling the partial wave have been suppressed. Demanding ^-T lowfc (/c', k; k 2 ) = implies 
an RG equation, which is claimed in Ref. and elsewhere to be 



d A _ 2 V A (k', A)T A (A, k; A 2 ) 

dK V {k ' k) -n 1 - (k/A) 2 • {2) 

Here we have switched to the abbreviated notation V A = Vowfc- By construction, the 
low-momentum HOS T-matrix of the initial interaction is preserved by the interaction V A 
obtained from integrating the RG equation. 

Reference [l| states that the RG equation, Eq. ([2]), must be modified in the event that 
bound states are present (i.e., in the 3 Si- 3 D! channel) to 
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The extra term arising from the bound states is given by [l[ 

where the bound state wave functions Xb-(^) = (^IXb ) satisfy 

(Ho + V A )\ X %) = -k 2 Bi \ X A Bi ), (5) 

and we have used the plane wave completeness relation P — - j A p 2 dp\p){p\ to obtain the 
second line in Eq. (J3J). 

We now provide a simple demonstration that the Sj3 A (k',k) term is identically zero. 
As a consequence of the cutoff independence of the HOS T-matrix, the interacting scat- 
tering eigenstates of the low-momentum Hamiltonian H A = H + V A are equal to the 
low-momentum projections of the corresponding eigenstates of the "bare" Hamiltonian, 
\Xk) = P\^k), with an analogous relation for bound states, \Xb) = P\^Bi) or equiva- 
lently Si(k'\V \Xb) — 0. Moreover, the preservation of the low- momentum part of the 
bare HOS T-matrix also implies that bound state poles of the input Hamiltonian should be 
preserved by H A (provided that the binding momentum is less than the sharp cutoff). With 
the aid of these observations we find 



d 

dA kBi 



-k% = , (6) 
d . a . 2 .„ _ , 2 



^Ixi) = -A 2 |A)(A|^) = -A 2 \A)(A\xi) ■ (7) 

Taking d/dA of both sides of the bound state Schrodinger equation and making use of the 
previous two equations leads to 

dV A 9 

-j^lxi) = ~-A 2 (k% +H + V A )\A)(A\ X %) . (8) 

Multiplying from the left by a P-space plane wave state (k'\ gives an expression for the first 
term in the brackets in Eq. (jlj), 

dV A 9 

(k'\^\xi) = --A 2 (k'\V A \A)(A\x%). (9) 

Finally, making the trivial substitution V A = H A — H in the second term in Eq. (j3J) gives 

_ 2 (k'\V A \A)(A\V A \x A m ) 2 (k'\V A \A)(A\(H A -H )\ X A Bi ) 

vr l + (k B JAy ' vr l + (k B JAy [ } 

9 :A 2 (A;'|V A |A>(A|;4> , (11) 
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which evidently cancels the first term. Therefore, we have shown that 5/3 A (k', k) = 0, and 
the original RG equation, Eq. (j2J), holds independent of the presence of bound states. The 
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equivalence of the model-space methods to the RG equation follows immediately, since it is 
well-known that such methods preserve the HOS T-matrix, the bound state poles, and the 
P-space projection of low energy eigenstates. (In addition, it follows directly that the extra 
terms in Eq. (4.8) of Ref. [IJ are cutoff independent.) By the same line of reasoning the 
smooth-cutoff version of the RG equation from Ref. [4j is also unmodified by bound states. 

Given the subtleties that can arise in any such demonstrations, a practical test is whether 
actual numerical applications of the RG equations reveal any problems related to complete- 
ness in channels with bound states. This test is not ideal because the numerical implemen- 
tation of the RG equations is difficult, and errors for phases shifts or bound-state energies 
grow in each channel as the cutoff is lowered. However, there is no sign of special problems 
or enhanced errors in the 3 Si- 3 Di channel [||. 

There are several other points in the manuscript that we feel need clarification or correc- 
tion. The first is the relationship of the RG evolution and chiral effective field theory (EFT) 
interactions. If one starts from a chiral EFT interaction determined to a given order in 
the power counting, the truncation error of the initial chiral EFT is exactly preserved when 
many-body interactions are evolved as well. If many-body interactions are truncated, we 
find that the error remains of the same natural size for cutoffs used in practical calculations. 
Therefore, V[ w/c preserves the systematic nature of the initial EFT. On the other hand, if 
one starts from a phenomenological interaction, then it is correct that the low-momentum 
Hamiltonian is not systematically improvable since the starting point is not either. 

There are also comments on the smooth-cutoff version of V\ ow k that may leave a wrong 
impression. In particular, the three-step numerical treatment applied in Ref. pj], which 
is used in practice, is numerically robust. While analytic analysis may be difficult, the 
converse is true for energy-dependent Wilsonian RG implementations, whose simplicity is 
highlighted in Ref. [l|. That is, while it is much easier to perform a fixed-point analysis 
within the energy-dependent framework, it is notoriously difficult to work with in few- and 
many-body calculations (to our knowledge no application to A > 3 exists). 

Finally, there are comments in Ref. |l| on the Similarity Renormalization Group (SRG) 
approach (inexplicably called WGW in Ref. [l[) about the non-invariance of the HOS T- 
matrix under the SRG, which seem to imply this poses a problem. We agree that the HOS 
T-matrix changes with SRG evolution, but disagree that this is a problem. The HOS T- 
matrix is certainly not an observable, while one can easily show that the on-shell T-matrix 
and bound-state poles are invariant under the SRG evolution. Indeed, the invariance of 
observables follows trivially from the fact that the SRG evolution represents a series of 
unitary transformations. 
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